Electronic viscosity in a quantum well: A test for the local density approximation 
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In the local density approximation (LDA) for electronic time-dependent current-density functional 
theory (TDCDFT) many-body effects are described in terms of the visco-elastic constants of the 
homogeneous three-dimensional electron gas. In this paper we critically examine the applicability 
of the three-dimensional LDA to the calculation of the viscous damping of 1-dimensional collective 
oscillations of angular frequency a; in a quasi 2-dimensional quantum well. We calculate the effective 
viscosity C,{u}) from perturbation theory in the screened Coulomb interaction and compare it with 
the commonly used three-dimensional LDA viscosity Y[uo). Significant differences are found. At low 
frequency Y{u)) is dominated by a shear term, which is absent in ^(tj). At high frequency C,{uj) and 
Y(uj) exhibit different power law behaviors and uj'^^"^ respectively), reflecting different spectral 
densities of electron-hole excitations in two and three dimensions. These findings demonstrate the 
need for better approximations for the exchange-correlation stress tensor in specific systems where 
the use of the three-dimensional functionals may lead to unphysical results. 



I. INTRODUCTION 

Time-dependent density functional theory is one of the 
premier techniques for the study of the dynamics of quan- 
tum many-body systems fii^ii^ The central idea of the the- 
ory, dating back to the pioneering work of Hohenberg, 
Kohn and Sham,**'^ is that the interacting many-body 
system can be simulated by a non-interacting system 
that yields the same density under the action of the self- 
consistent electrostatic potential Vu and an additional 
exchange-correlation (xc) potential Vxc that is uniquely 
determined by the particle density This approach, 
while rigorous in principle, runs into serious difficulties 
when applied to phenomena that involve dissipation and 
memory - for example, the damping of collective modes. 
This is because Vxc is a strongly non-local functional of 
the density!^ 

About 20 years ago a time-dependent current den- 
sity functional theory (TDCDFT) was proposed)^^ which 
offered a natural way to treat memory and dissipa- 
tion without losing the advantages of a local descrip- 
tion of many-body effects. In this theory the exchange- 
correlation potential was replaced by an exchange- 
correlation force field F^c = —eExc (the notation is de- 
signed for electronic systems: — e is the electron charge 
and Exc is an "exchange-correlation electric field" ) which 
was represented as the divergence of a stress tensor a 
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which in turn was expressed as a local functional of the 
equilibrium density n and the velocity u = ^, where j is 
the particle current density. 

In the linear response regime an approximate form of 
the stress tensor was proposed which has the same 
form as the stress tensor of the classical Navier-Stokes 
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where pxc is the exchange-correlation pressure of the ho- 
mogeneous electron gas at the local density n, rj is the 
shear viscosity and ^ is the bulk viscosity of the homoge- 
neous electron gas at the local equilibrium density uq. rj 
and C, are frequency-dependent and have imaginary parts 
which are related to a shear modulus and a dynamical 
bulk modulus respectively. 

An essential feature of this approximation is that it is 
still local in the equilibrium density and the velocity field. 
In particular, the adiabatic-local density approximation 
(3D-LDA) amounts to keeping only the first term on the 
rhs of Eq. dH). 

The exchange-correlation field ([T|) with stress ten- 
sor has been applied to the study of sev- 
eral systems (semiconductor quantum wells j^^'^^i^^ 
atomSfii semiconductors and polymerS j^^'^^i^^ molecular 
junctionsiSiiSiii^iSS and metal o^^'^^ ) with varying degrees 
of success. However a fundamental difficulty exists. Be- 
cause the viscosities rj and C are borrowed from an in- 
finite three-dimensional electron gas, this approximation 
implicitly assumes the existence of a continuum spectrum 
of excitations (electron- hole pairs). Needless to say, this 
assumptions is not justified in systems with discrete en- 
ergy levels. This may lead in some cases to spurious 
results. For example, the optical transitions between dis- 
crete energy levels of atoms are found to have a spurious 
linewidth.^^« 

A more "friendly" system is the electron gas in a two- 
dimensional quantum well (shown in Fig. [T]), in which 
the electrons perform collective oscillations in a direction 
perpendicular to the plane of the quantum well (the z- 
direction), while remaining homogeneous in the plane of 



2 



QW 



2DEG 



v{r) 

-< ► 

Coulomb 
interaction 



- "v / 



1. ■ 



FIG. 1: Two-subsystems model for electrons in a two- 
dimensional quantum well. The collective motion of the elec- 
trons in the z direction is coupled, through the Coulomb in- 
teraction, to density fluctuations in the plane of the 2DEG. 
Energy is exchanged between the two subsystems via the cre- 
ation of electron-hole pairs in the 2DEG. 



the well (the x — y plane). Spontaneous oscillations of 
this kind - known as intersubband plasmons - have been 
intensively studied because of their possible use in de- 
vices in the terahertz frequency range. In practice these 
oscillations, spontaneous or forced, provide one of the 
best testing grounds for TDDFT calculations of damp- 
ing and relaxationJ^iS^i^i Even in this favorable case. 



however, one may question the accuracy of the local den- 
sity approximation for the viscosity. To begin with, the 
electron-hole pairs that are responsible for the damp- 
ing and the screening of the electron-electron interaction 
should be those of an essentially two-dimensional elec- 
tron gas (2DEG) - not those of a three-dimensional elec- 
tron gas as implicitly assumed in the 3D-LDA. However, 
the system is not exactly two-dimensional, even if one 
neglects inter-subband transitions, because the electron- 
electron interaction in the lowest subband is modified by 
form factors, which take into account the finite extent 
of the wave function in the z direction. More impor- 
tantly, in a collective oscillation that preserves uniformity 
in the x ~ y plane, the Fermi surface of the 2DEG does 
not change its (circular) shape (see Fig. [2]): this means 
that only the hulk viscosity, contributes to the damp- 
ing of the collective mode. By contrast the dominant 
contribution to the damping in 3D-LDA comes from the 
shear viscosity, which is associated with changes in shape 
of the implicitly assumed three-dimensional Fermi sur- 
face. Finally, it is evident that the above difficulties can- 
not be solved by resorting to a strictly two-dimensional 
LDA, since the system becomes homogeneous and time- 
independent when strictly projected in the x — y plane. 
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FIG. 2: (a) In a three-dimensional collective mode the Fermi surface changes shape periodically from oblate to prolate. 
The process can be pictured as a transfer of electrons between the Fermi disks into which the Fermi surface can be sliced, 
perpendicular to the direction of propagation of the mode. These excitations are responsible for the finite shear viscosity of 
the electron liquid, (b) In 2D, for an isotropic single band system the Fermi surface is a circle: excitations can only change the 
radius of the circle, not its shape. There is no shear viscosity, (c) In a multi-band isotropic 2D systems, inter-band transitions 
allow a change in shape of the Fermi surface - formed in this case by several parallel disks: the shear viscosity reappears. (2D 
wave vectors are represented by capital letters). 



In view of the above difficulties it is clearly of great in- 
terest to test the validity of the 3D-LDA against a more 
direct calculation of the damping rate of a collective os- 
cillation - a calculation that explicitly takes into account 
the reduced dimensionality of the quantum well. Our 
idea is to study a particularly simple oscillation that is 
described, to zero-order in the Coulomb interaction, by a 



separable wave function, i.e. the product of an assigned 
time-dependent function of the z coordinate of the elec- 
trons times a function of the x and y coordinates: the 
latter may be either the ground state or a uniform ex- 
cited state of the quasi-2DEG in the x — y plane. For 
example, we can put all the electrons in a given time- 
dependent combination of the lowest and the first ex- 
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cited subband: we assume that there is an external time- 
dependent potential that can do this. The main physical 
assumption here is that all the electrons populate a single 
time-dependent subband, and no inter-subband scatter- 
ing is considered. The time-dependent electronic density 
n{z,t) creates, via time-dependent form factors, a time- 
dependent Coulomb field which induces transitions be- 
tween different states of the 2DEG. We can make use of 
the Fermi Golden Rule to calculate the rate at which the 
motion in the z-direction creates electron-hole pairs in 
the x^y plane. From this we calculate the rate at which 
energy is transferred from the oscillatory motion in the 
z direction to electron-hole pair excitations in the x — y 
plane. In brief, the electron-hole pairs in the x — y plane 
act as a thermal bath for the motion in the z direction. 
Obviously this approach is justified in the limit of weak 
Coulomb interaction (high density). 

It turns out that it is possible to recast the results 
for the damping rate in a form that is analogous to the 
LDA, except that it involves a non-local viscosity z'). 
This non-local viscosity can be converted into an effective 
local viscosity C,{z) by integrating over z' at fixed z - a 
procedure that amounts to neglecting the finite range 
of nonlocality on the scale of variation of the density. 
In this manner we are able to compare the approximate 
LDA viscosity Y{z) - a function of the local equilibrium 
density - with the more accurate one, C,{z). 

We find significant differences between (^{z^uj) and 
Y{z,uj). At low frequency the latter is dominated by a 
shear term, which is absent in C,. For this reason C, turns 
out to be numerically much smaller than Y . At high fre- 
quency C,{z,Lij) and Y{z,ijj) exhibit different power law 
behaviors (w"'^ and uj~^l'^ respectively), reflecting differ- 
ent spectral densities of electron-hole excitations in two 
and three dimensions. 

We evaluate the damping rate as a function of the 
frequency of the external field, w, for the screened and 
unscreened potentials. We find, as one could expect, 
that the rate of energy dissipation in the screened (un- 
screened) interaction goes as (w^) at small frequen- 
cies. In the 3D-LDA approximation, instead, the energy 
dissipation rate goes as cj^. Thus the 3D-LDA largely 
overestimate the energy dissipation. 

It is already knowrJ^ that the LDA viscosity is spu- 
rious in small systems like atoms and molecules. The 
present findings show that a serious loss of accuracy can 
occur also in infinite systems, when the motion is re- 
stricted to a single one-dimensional subband as in the 
example considered here. Therefore, one needs to ap- 
ply special care when using 3D-LDA approximations in 
specific systems such as the one we consider here. The 
development of a better approximation for the viscous 
stress tensor - an approximation that is uniformly appli- 
cable across different dimensionalities - thus emerges as 
a critical issue. 



II. THE MODEL 

The Hamiltonian for an interacting N-electron system 
confined to a quantum well (QW) of width L in the z 
direction is 

where v[r) — e^/47re|r| is the Coulomb potential (e is the 
background dielectric constant) and 

In the absence of electron-electron interactions (y(r) = 
0) we can separate the Hamiltonian into two parts de- 
scribing respectively the motion in the z direction and 
the motion of the two-dimensional electron gas (2DEG) 
in the plane: 

Ho^H,+ H2DEG (5) 

where 

N 2 

^- = E§J + ^(^») + ^-*(^'*)' (6) 

i=l 

N 2 , 2 N 2 

H2DEG-l^ 2m =2^2^^ 

i—l i—1 

where Vext{z,t) is the external potential that excites 
the QW (see the following discussion). Obviously the 
Coulomb interaction couples these two subsystems in a 
non-trivial way. In the following we will consider the 
Coulomb interaction as a small perturbation and then 
apply time-dependent perturbation theory to study the 
energy exchange between the z-motion and the density 
fiuctuations of the 2DEG. In particular, we will show 
that, if one assumes a given motion in the z direction, 
the Coulomb interaction excites electron-hole pairs in the 
X — y plane, thus effectively transferring energy from the 
z-motion to the 2DEG. 

A solution to the Schrodinger equation with the unper- 
turbed Hamiltonian Hq can be written as the product of 
solutions of H2DEG and in the following manner: 

I*) = \^)2Deg\x). (8) 

where \^)2DEG is a Slater determinant of N plane waves 
in the coordinates Ri and 

\x)z - X |x)2... X \x)n (9) 

is a symmetric function of the z coordinates of the elec- 
trons. We assume that all the electrons share a common 
motion in the z direction specified by the one-electron 
wave function x(z, <), which is sustained by a suitable 
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time-dependent external potential Vext{z,t). The den- 
sity remains uniform in the x — y plane. The antisym- 
metry of the complete wave function is ensured by the 
antisymmetry of |$)2_d_bg- 

Next we turn on the electron-electron interaction, 
and derive an effective hamiltonian for the dynamics 
of \^)2DEG in the presence of the driving potential 
Vext{z,t). This is done by substituting Eqs. ^ and 
^ into the Schodinger equation 



ihdt\^{t)) =H{t)\^{t)). 



(10) 



Making use of the Fourier representation of the Coulomb 
potential, 



q,Q 
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{R and Q are vectors in the 2D real and momentum space 
respectively, with r = {R,z), g is a ID momentum and 
V the volume), and the Fourier transform of the density 
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we can write the electron-electron interaction as 



i=^i q,Q 



(12) 
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Assuming that the total wavefunction is still approxi- 
mately of the form ([5]) we can ignore the correlation be- 
tween the z's and the i?'s, and simply replace e"**'^'"^^-* 
by its average value in the given one electron wave func- 
tion dH): 
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= \n{q,t)\^ 



(14) 



independent of i and j. Then the effective dynamics for 
\^)2DEG is given by 



ihdt\^)2DEG = Heff\^)2DEG 



(15) 



where 
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Notice that |n(q,t)p plays the role of a time- dependent 
form factor for the electron-electron interaction in the 
plane. 



To apply perturbation theory we assume that 7i(g, t) is the sum of a static part plus a small dynamical contribution: 

n{q,t) = no{q) + Sn{q,t) (17) 

with \Sn{q,t)\ <^ ?io('z) for all q and t. Neglecting contributions of order ((5n)^ we get the effective time-dependent 
Schrodingcr equation for |$)2_d_bg 
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The left hand side of Eq. (fT8|) is the sum of a static part 
and a time-dependent part. The static part includes the 
effect of the QW static form factor on the ground state 
and the excited states of the 2DEG. We assume that 
the 2DEG is in its ground state at t = 0. The time- 
dependent part induces transitions between eigenstates 
of the static 2DEG Hamiltonian and thus is responsible 
for dissipation. 

An important check for this model is that if the mo- 
tion in the quantum well is a rigid translation of the 
equilibrium density in the z direction (Kohn mode), i.e., 
n{z, t) = no{z + u{t)) ~ 71,0(2;) +dzno{z)u(t) where u{t) is 
a spatially uniform function of time, then the dissipation 
must be absent. This is easily verified from Eq. 
for the Kohn mode we have Sn{q,t) = iqno{q)u{t) where 



u{t) is a real function, and 

Re[no{q)STi{-q,t)] = Re[iq\no{q)\^u{t)] = . (19) 

Therefore only the static non-dissipative part of the 
Hamiltonian survives for the Kohn mode. 

III. PERTURBATION THEORY 

To apply the Fermi golden rule we assume that the 
forcing density 5n(z, t) is a periodic function of time with 
angular frequency oj, Sn{q,t) = 6n{q,uj)e~'"^* + c.c, and 
we write 

= (20) 
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where j{q,uj) is the Fourier amplitude of the current 

no{z)dtu{z, t) = j{q, c^)e^(?^-*) + c.c. (21) 

The energy transferred per unit of time by the 2:-motion 
to the 2DEG is equal to the energy huj absorbed in each 
allowed transition from the ground-state to an excited 
state of the 2DEG times the rate at which the transition 
occurs. Evidently the only allowed transitions in lead- 
ing order are the ones in which two electron-hole pairs 



with opposite wave vectors Q and —Q are created in 
the 2DEG, with energies adding up to huj. The first 
pair is created by promoting an electron from state K 
within the two-dimensional Fermi surface of the 2DEG 
to state K + Q outside the Fermi surface. Similarly, the 
second pair is created by promoting an electron from K' 
to K' — Q, with K' inside and K' — Q outside the Fermi 
surface [see Fig. [5J panel b)] . So the energy absorbed per 
unit time is 
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(22) 



where tik is the Fermi occupation number [tik = Q{Kf — K), where G(x) is the step function and Kp is the Fermi 
wave vector] and ujk,k',q = fi[Q^ + Q ■ {K' — K)]/m is the energy of the double electron-hole pair. A is the area 
occupied by the 2DEG. The two terms in the square bracket arise from direct and exchange processes respectively - 
the latter arising from the antisymmetry of the 2DEG wave function. In the following we will disregard the exchange 
term. Aside from the fact that its contribution is small in the high density limit, we must keep in mind that the 
exchange contribution was also dropped in previous calculations of the LDA viscosity: ■^^ for a fair comparison we must 
drop it here too. 

Introducing in Eq. (|22p the imaginary part of the two-dimensional Lindhard function- 
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(23) 



where ojk.q = ft[Q^/2m — Q ■ Kim] we obtain 
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(g2+Q2)(^/2+Q2) 



Notice that Xo^ depends only on the sheet density of the electron gas, which is independent of z. 
As a last step we express this result in terms of the velocity field, v[z, t) = dtu{z, t) = j(z, t)/no(z). 



(24) 
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where 



C(z,z';w) =no(z)no(z')^ / d^' I dQ Qv{Q, z)i* [Q, z')'^mxl^ {Q,oj')'^mxl° {Q,oj - u') 



(26) 



is the non-local bulk viscosity of the electron gas in the QW and 



v{Q,z) 
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dq 
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(72 + Q2 



(27) 



the effective electron-electron interaction in the plane of the 2DEG. 
Eq. (|25p should be compared with the corresponding LDA expression 
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(28) 



6 



where (see Refs. |25| and [2; 
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Xo^(9,w) is the 3D Lindhard function, Xoril^^) noninteracting transverse current-current response function 

in 3D, and v{q) is the Fourier transform of the screened 3D Coulomb interaction (about the screening more will be 
said later). The z-dependence of these functions arises from the fact that the 3D Lindhard functions depend on the 
local density. 

To facilitate the comparison between the two sets of formulas p5 |) - ([26)) and ([28l) - ((30)) we introduce the integrated 
viscosity 

az,u)= f C{z,z';co)dz' , (31) 



which has the same physical dimensions as (^^^^{z,uj) [Energy x Time/Volume] and can be directly compared to it. 
The explicit expression for C(z,a;) is 

C{z,uj)=noiz)^ duj' J dQ Qi{Q,z)v* {Q)Qmxl''{Q, uj')^fnxl''{Q, cu-iv') , (32) 

where 

v{Q) = / no{z)v{Q, z)dz . 



It is interesting to compare Eq. with the expres- 

sions obtained from the xc kernel in the 3D-LDA [see 
Eq. ([29|l]. In the homogeneous 3D electron liquid there 
are two independent viscosities: the bulk viscosity C([j) 
and the shear viscosity r]{uj). The former arises from mo- 
tions that change the volume of the local Fermi surface 
(the density), but not its shape. The latter appears when 
the motion changes the shape of the Fermi surface, even 
if its volume does not change. A coUisionless longitudinal 
collective mode, such as the ordinary plasmon, involves 
both types of motion simultaneously, as discussed in the 
caption of Fig. [51 the effective viscosity for such a mode 
isY{w) = C^"^{cj) + j7f^^{u;)^^Ycanhe easily con- 
structed from the combination of the two equations (|29|) 
and ([30ll . 

Aside from the obvious difference in the dimensionality 
of the Lindhard functions, we see that our expression 
for ({z,uj) is formally similar to the expression ((29l) for 
C^^^^iuj) in 3D-LDA. These expressions, however, van- 
ish as at low frequency, because the Lindhard spec- 
tra Smxo^^^^^ ((/, w) vanish as w. As a result ({z,uj), 
as well as C^^^(a;), goes as at low frequency. On 
the other hand, the behavior of the LDA shear viscosity, 
ri^^^{uj), is quite different. From Eq. (|5D)) we see that 
this quantity contains a term involving the transverse 
current spectrum Smxp-r j ^^^^ term tends to a finite 



limit for cu ^ because the factor l/cu^ compensates for 
the smallness of the Lindhard spectra. 

Therefore a fundamental difference exists between our 
results and those of the 3D-LDA in the low frequency 
regime. While the LDA viscosity is dominated by the 
shear term, which remains finite for — s- 0, the present 
viscosity is purely of the bulk type and vanishes as 
for w ^ 0, at least when the screening of the electron- 
electron interaction is properly taken into account (see 
discussion below). Physically, the absence of a shear 
term is due to the fact that the oscillatory motion in 
the z direction does not change the shape of the local 
Fermi surface of the 2DEG (a circle). Whereas, in a 3D 
oscillation the local Fermi surface changes its shape pe- 
riodically from a prolate to an oblate ellipsoid (passing 
through the sphere), generating shear friction in the pro- 
cess (see Fig. ^ 

Before proceeding to a detailed comparison of the nu- 
merical results for the viscosities and the energy dissipa- 
tion rates we need to say something more about the role 
of the screening of the electron-electron interaction. First 
of all, it must be noticed that the integrals over wave vec- 
tor in Eqs. (glll-lISni) diverge, in the limit w ^ 0, if v{q) 
is taken to be the Fourier transform of the bare Coulomb 
interaction, Ane^/eq'^. The divergence comes from the 
small q region. The standard cure for this type of diver- 
gence is to screen the interaction by the static dielectric 
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function of the electron gas, which, in the high-density 
hmit, can be reasonably approximated by the Thomas- 
Fermi formula. The screened interaction has the form 



v{q) 



(33) 



where k = ^y4^Te^N(0) is the Thomas-Fermi wave vector, 
and N{Q) is the density of states at the Fermi energy. 

The situation is somewhat different in Eq. p2|) . Be- 
cause the effective interaction v{Q,z) diverges at small Q 
only as , it turns out that the wave vector integral is 
finite even without including the screening. Nevertheless, 
the integral has a large contribution from extremely small 
wave-vectors, of order uj/vp for w — > 0, which effectively 
changes the formally expected behavior of ({z,u!) to 
a constant independent of to. This curious phenomenon 
is shown in Fig. |4l where we see that ({z,uj) tends to a 
constant when the interaction is not screened. Inclusion 
of screening however, will drastically modify this behav- 
ior, reinstating the expected behavior (see Fig. [5]). To 
demonstrate this point we have calculated (^(z^uj) with 
the bare interaction v{Q, z) replaced by a screened inter- 
action according to the scheme 



v{Q,z) 



v{Q,z) 



i-v{Q)xo{Q,0) 



(34) 



The denominator is the static dielectric constant of a 2D 
QW in which all the electrons reside in the lowest sub- 



band, and no intersubband transitions are allowed. In- 
deed the low-energy excitations of the 2DEG are exclu- 
sively mira-subband electron-hole pairs. Therefore ne- 
glecting mter-subband transitions is justified at low fre- 
quency and qualitatively correct at higher frequencies (as 
long as only a few subbands are involved). 



IV. NUMERICAL RESULTS 

Fig. [3] shows the calculated values of the non-local vis- 
cosity C(z, z', (jj) as a function of z and z' for a sheet den- 
sity n2D = 10"'^'^ cm~^ and L — 40 nm {Kp — \J2-Kn2D — 
7.9 X 10''^ cm~^ ) with and without the screening (pan- 
els a) and b) respectively) . We have assumed that the 
equilibrium density in the QW has the form 



no(^) 



(?)■ 



appropriate to the lowest subband of a free particle in a 
box. This choice gives 



v{Q,z) = - 



27re2L / e"^ cosh(Qz) cos 



and, for the unscreened interaction. 



- 1 



(35) 
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(36) 



The functions v{Q^ z) and v{Q) depend on the physical parameter KpL where Kp is the Fermi momentum of the 
2DEG and L is the thickness of the well in the z direction. For KpL ^ 1, v{Q,z) ~ v{Q) ~ 2TTe^/eQ, i.e. the 
functions v reduce to Fourier transforms of the bare 2-dimcnsional Coulomb potential. 



The nonlocal viscosity ^(z, z', a;)/(rio(z)?i) decays when 
one of the points is close to the boundary and it is maxi- 
mum when both z and z' are near the center of the well. 
While the spatial behavior of ({z,z',uj) looks qualita- 
tively similar for the screened [panel a)] and unscreened 
case [panel b)], the order of magnitude is quite different, 
the unscreened results being quite large. This has its 
origin in the strong singularity of the Coulomb potential 
for small momenta. We then expect that in the case of 
the unscreened interaction, the dissipation rate be largely 
overestimated. 

In Fig. [3] we plot the integrated viscosity C{^^^) for 
various values of KpL as a function of the frequency for 



the unscreened interaction. For all the values of the pa- 
rameter KpL we notice the existence of a finite limiting 
value for small frequencies. As already discussed, this 
finite limiting value is due to the existence of a strong 
singularity for small Q, l/Q"^, which compensates for the 
smallness of the Lindhard spectra. Including the screen- 
ing according to the scheme of Eq. (|34p changes this be- 
havior dramatically. The new behavior is shown in Fig. [5] 
for various values of the parameter KpL. Now C,{0,uj) 
vanishes as in agreement with the general discussion 
of the previous section. 

Finally in Fig. [S]we compare the behavior of C.{0,uj) 
with the corresponding LDA quantity Y{u}), evaluated at 
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FIG. 3: The nonlocal viscosity (^(z, z' ,u!)/no as a function of z, z' for different values of with KpL — 3.17. Panel a) - screened 
interaction; panel b) unscreened interaction. Notice the different scales in the two panels. 
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FIG. 4: Average local viscosity (^{z,uj)/hno{z) from Eq. (|32p 
as a function of frequency for z = and for different values 
of KpL. The Coulomb interaction potential is unscreened. 
Notice that approaches a finite value for cj —> 0. 



the density nQ{z = 0). In the screened case our viscosity 
falls well below Y{lu) (calculated from the parametriza- 
tion of Qian and VignaleS^) demonstrating the intrinsic 
limitation of the 3D-LDA in this regime. 

It is also interesting - although not fully justified phys- 
ically - to consider the high-frequency behavior of C{z, 
vis-a-vis the high-frequency behavior of Y(lu). This can 
be calculated analytically and one sees that the relevant 
contribution to the integral over wave vector comes from 
Q proportional to V 2mw / h and is therefore large. In 
this regime both v{Q) and v{q) scale as uj~^, so the only 
difference between the expressions and (|29p comes 
from the different "volume element" in wavevector-space, 
Q and respectively. Taking this into account we im- 
mediately understand the origin of the ^/lu difference in 
the high-frequency behaviors of (^(z,uj) {ll>^^) and Y{ll!) 
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FIG. 5: Same as Fig. U but now for the case of a screened 
Coulomb interaction. Notice that ("(0, lj) oc for to ^ (the 
dashed line is the slope of lj^ in the logarithmic plot). 



(a;~^/^) - the former tending to zero faster than the lat- 
ter. 

To fully evaluate the energy dissipation via either 
Eq. p5|) or Eq. ((28)) we need to assign the velocity v{z, uj). 
In our model we consider coherent oscillations of the wave 
function in the quantum well. For this motion we then 
choose 



(37) 



where A is a constant and Xn{z) are the normalized eigen- 
functions of the quantum well 



Xniz) 



L 



n even 



L 1 cos (2^) n odd 



(38) 



For the linear response theory to be valid we have as- 
sumed A ^ li^ We can evaluate the current density for 
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the state and obtain the velocity field 



d;:V[Z,Uj) = dz -T-- = -UJ\X\ COS — — 



(39) 



FIG. 6: Comparison between (^{z,ij)/hno{z) and the 3D-LDA 
viscosity Y{z,uj) for z — 0. As expected the screened and 
unscreneed results for ({z,uj) coincide for large frequencies. 
Notice that the 3D-LDA (with a screened interaction) predicts 
a finite limit of the viscosity for lj ^ 0. 



By substituting this expression in the equation for the 
energy dissipation, Eq. ((25|) we get 



dE _ Ah 

dt TTUJ 



dio' / dQ Q3mxr(Q,w')Smxr(Q-^-^')r'(g,w') 



(40) 



where we have defined for the unscreened case 



r(Q,c^)-- 



e^uj \X\L 
27re 1 . 



e 



¥ sinh 



1 - 



[LQY 



2QL 



1 (-^9)" 



(41) 



Accordingly to the scheme we adopted [see Eq. ([M)) ]. to 
include screening in the interaction potential we substi- 
tute r{Q,uj) with T{Q,oj)/[l - Xo{Q)v{Q)]. 

In Fig. El we plot d£/dt = d/dt{E/Nhiu^\X\^), where 
ujp — 2Trno{Q)e^ Kp/m is the 2D plasmon frequency eval- 
uated at the Fermi momentum, by directly evaluating 
Eq. (PO]) for the screened and unscreened case. It is 
immediately seen that for small frequencies the energy 
dissipation rate scales as (w^) for the screened (un- 
screened) interaction. This has to be compared with the 
energy dissipation rate obtained from the 3D-LDA: since 
Y{lu) goes to a constant for small frequencies [see Fig.[B], 
we see that the 3D-LDA dissipation rate scales as oj^ for 
small oj. Then we can conclude that the 3D-LDA over- 
estimates the energy dissipation rate in our model. 

V. DISCUSSION 

The system studied in this paper can be thought of as 
consisting of two weakly coupled sub-systems of reduced 
dimensionality: a linear oscillator in the z direction cou- 
pled to a two-dimensional electron gas in the x — y plane. 
We have calculated the energy transfer between these 



two sub-systems to leading order in the strength of their 
coupling and we have thus identified a non-local viscos- 
ity which can be compared (after the integration of one 
variable) with the viscosity obtained from the standard 
3D-LDA. We have found very significant differences be- 
tween the two viscosities, particularly at low frequency, 
where the 3D-LDA viscosity is dominated by a shear term 
which is absent in the present treatment. 

It is important to realize that the problem we have 
uncovered stems from the extreme one-dimensional char- 
acter of the dynamics of our model. We have assumed 
that all the electrons reside in a single time-dependent 
subband x(z,t). In effect, this dramatically limits the 
type of excitations we can generate in the 2DEG (see 
Fig. [2]) and this is the reason why we do not have the 
transverse current term that is responsible, in three di- 
mensions, for the shear viscosity term. It is expected - 
and this case is currently under investigation - that a 
transversal shear viscosity term will reappear in a multi- 
subband system via the relative motion of the different 
sheets of the Fermi surface in different subbands. For a 
strictly one-dimensional motion, however, it seems clear 
that the use of the LDA is problematic. This leaves us 
with the challenge of formulating a better approximation 
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to keep track of the systematic reduction in viscosity fol- 
lowing from the geometric confinement of the system. 



FIG. 7: The energy dissipation rate for the screened and un- 
screened (inset) interaction [see Eq. (|40|l ] as a function of the 
external frequency for various values of the quantity KpL. 
The energy dissipation rate scales as oj'* for the screened inter- 
action, a result much smaller than the uj'^ behavior obtained 
from the 3D-LDA energy dissipation rate (also for a screened 
interaction). The limiting behavior d£/dt ~ a;* is shown for 
reference as a dotted line. Inset: Due to the strong singular- 
ity of the unscreened interaction, the energy dissipation rate 
for this case, scales as uj'^ . The dotted line shows the limiting 
behavior d£ / dt ui^ . 
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